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Abstract 

We calculate the 0(a s ) radiative corrections to polarized top quark decay into a 
charged Higgs boson and a massive bottom quark in two variants of the two-Higgs- 
doublet model. The radiative corrections to the polarization asymmetry of the decay 
may become as large as 25%. We provide analytical formulae for the unpolarized 
and polarized rates for nib ^ 0 and for mb = 0. For mb = 0 our closed-form 
expressions for the unpolarized and polarized rates become rather compact. 


1 Introduction 


The purpose of this paper is the evaluation of the first order QCD corrections to the decay 
of a polarized top quark into a charged Higgs boson and a bottom quark. Highly polarized 
top quarks will become available at hadron colliders through single top production, which 
occurs at the 33 % level of the top quark pair production rate [HEj. Future e + e _ - colliders 
will also be copious sources of polarized top quark pairs. The polarization of these can 
easily be tuned through the availability of polarized beams Pi]. Measurements of the 
decay rate and the polarization asymmetry in the decay t(j) —> H + + b will be important 
for future tests of the Higgs coupling in the minimal supersymmetric standard model 
(MSSM). 

The 0(a s ) corrections to the decay rate t —> H + + b have been calculated previously 
in 0 [6j, [H (8, [9, TU| [TT} [U, 13], and have been found to be important. The present 
paper provides the first calculation of the 0(a s ) radiative corrections to the polarization 
asymmetry in polarized top quark decay f (j) —> H + + b. Depending on the value of tan /3 
and the mass of the charged Higgs boson the radiative corrections to the polarization 
asymmetry can become quite large (up to 25 %) in one variant of the two-Higgs-doublet 
model and must therefore be included in a decay analysis. 

The decay t —> H + + 6 is analyzed in the rest frame of the top quark (see Fig. 1). The 
three-momentum q of the H + boson points into the direction of the positive z-axis. The 
polar angle Op is defined as the angle between the polarization vector P of the top quark 
and the z-axis. 



Figure 1: Definition of the polar angle dp. P is the polarization vector of the top quark. 


lit terms of the unpolarized rate T and the polarized rate T p the differential decay 
rate is given by 


^ = i (r + r r- co S e P ) = \ r (i + cos e P ). (U 

where P is the degree of polarization. The polarization asymmetry ap is defined by 
ap = r p /r. Alternatively one can define a forward-backward asymmetry A FPj by writing 
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( 2 ) 



where F F and T F are the rates into the forward and backward hemispheres, respectively. 
The two measures are related by A F b = \Yau- In our numerical results we shall always 
set P = 1 for simplicity. 

Technical details of our calculation can be found in [14]. As in [14] we use dimensional 
regularization (D = 4 — 2u with w <C 1) to regularize the ultraviolet divergences of the 
virtual corrections. We regularize the infrared divergences in the virtual one-loop and the 
real gluon emissions (tree-graph) corrections by introducing a finite (small) gluon mass 
m g 0 in the gluon propagator. In the tree-graph corrections the phase space boundary 
becomes deformed away from the IR singular point through the introduction of a (small) 
gluon mass. The logarithmic gluon mass dependence resulting from the regularization 
procedure cancels out when adding the virtual and tree-graph contributions. We have 
checked consistency with the Goldstone boson equivalence theorem which, in the limit 
mw+ /m t —■► 0 and m F + /m t —> 0, relates Y and Y p for t —> H + + b to the unpolarized and 
polarized longitudinal rates Yl and Y p in the decay t —> W + + b calculated in [14] . 

2 Born term results 

The coupling of the charged Higgs boson to the top and bottom quark in the MSSM can 
either be expressed as a superposition of scalar and pseudoscalar coupling factors or as a 
superposition of right- and left-chiral coupling factors. The Born term amplitude is thus 
given by 



If the top quark mass rn t , bottom quark mass mb and the charged Higgs boson mass 
m H + satisfy m t > m H + + mb, the top quark decay t —> H + + b is possible in the MSSM. 


In order to avoid flavor changing neutral currents (FCNC) the generic Higgs coupling to 
all quarks has to be restricted. In the notation of [F5] in model 1 the doublet H i couples 
to all bosons and the doublet H 2 couples to all quarks. This leads to the coupling factors 
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model 1: 


a = 


9w 

2\f2rn w 


Vtb(m t - m b ) cot (3, 


b 


9w 

2y/2rn w 


V t b(m t + m b ) cot (3, 


(4a) 

(4b) 


where Vtb is the CKM-matrix element and tan (3 = Vi/v\ is the ratio of the vacuum 
expectation values of the two electrically neutral components of the two Higgs doublets. 
The weak coupling factor g w is related to the usual Fermi coupling constant G F by g 2 = 

\\[2m 2 w G F . 

In model 2, the doublet Hi couples to the right-chiral down-type quarks and the 
doublet H 2 couples to the right-chiral up-type quarks. Model 2 leads to the coupling 
factors 


model 2: 


—— Vt b (m t cot (3 + m b tan (3), 

2\J2rn w 

(5a) 

—^— V tb (m t cot (3 — rn b tan (3). 

2 V2m w V 

(5b) 


Since rn b <C rri t the b mass can always be safely neglected in model 1. In model 2 the 
left-chiral coupling term proportional to rn b tan [3 can become comparable to the right- 
chiral coupling term m t cot f3 when tan (3 becomes large. One cannot therefore naively 
set m b = 0 in all expressions in model 2. One has to distinguish between the cases in 
which the scale of m b is set by rn t cot 2 (3 and those in which the scale of rn b is set by 
mt . In the former case one has to keep rn b finite. In the latter case m b can safely be set 
to zero except for logarithmic terms proportional to In (m b /mt) that appear in the NLO 
calculation. Keeping this distinction results in very compact closed-form expressions for 
the radiatively corrected nnpolarized and polarized rates which we shall list in the main 
text. We shall, however, also present general m b ^ 0 results along with the rn b = 0 results. 
It is quite evident that one has to use the m b ^ 0 expressions for charged Higgs masses 
close to the top quark mass. 

For the amplitude squared one obtains 

\M 0 \ 2 = ^2 MlMo = 2(p t • p b )(a 2 + b 2 ) + 2 (a 2 - b 2 )m t m b + 4 abm t (p b ■ s t ), (6) 

Sb 
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where p t and pb are the four-momenta of the top and bottom quarks. s t denotes the 
polarization four-vector of the top quark. Accordingly the unpolarized and polarized 
Born term rates read 


r 


Born 


r p 

Born 


m t V\ 

167T 


(a 2 + b 2 )(l — y 2 + e 2 ) + 2(a 2 



m t a/A 
167t 


2abV \, 


(7a) 

(7b) 


where we have used the abbreviations A = A(l, y 2 , e 2 ) with the Kallen function A(a, b, c ) : = 
a 2 + b 2 + c 2 — 2(a6 + be + ca). The scaled masses y are e are defined as 


rn,H+ mb 

V :=-, e:=—• 

m t m t 


( 8 ) 


As mass values we take rrib = 4.8 GeV and m t = 175 GeV such that e = 0.027. For later 
use we register that the factor PS 2 = VX/(16nmt) is a two-body phase space factor. 

In the limit mb —> 0 the unpolarized and polarized Born term rates simplify to 


lim r Born 

m;,-* 0 


, o , 2Ul a 2 — b 2 2e - 

(a + b ){l + 2 2 --^}r, 

a z + b z 1 — y z 1 


lim T 

rrib^O 


Born 


= 2abT, 


where 


2\2 


r = 


m t ( 1 - y 2 ) 
167T 


(9a) 

(9b) 

( 10 ) 


Eqs. filial) and (l9bl) are quite useful for a discussion of the qualitative behavior of the 
unpolarized and polarized rates in the two models as long as one is not too close to the 
kinematical limit y — 1 — e. In fact, for most of the y region away from the endpoint 
y = 1 — e, Eqs. (j9aj) and (l9bl) are an excellent approximation to Eqs. (17a|) and (l7bl) . 

For model 1 one has (a 2 — b 2 )/(a 2 + b 2 ) = —2e/(l + e 2 ) and thus the second term 
in the curly brackets of Eq. fl9al) can be safely dropped for e —> 0 except when y gets 
close to 1 — e. For model 2 one has ( a 2 — b 2 )/(a 2 + b 2 ) = 2e/(cot 2 (3 + e 2 tan 2 (3), which, 
when combined with the factor 2e in Eq. f[9a]) . can become as large as 0(5%) for e.g. 
tan (3 ~ 5. Since this is further multiplied by (1 — y 2 ) -1 we shall therefore always keep the 
second term in the curly brackets of Eq. (l9a|) in model 2 when making use of the mb —■► 0 
approximation. 
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Fig. 2a Fig. 2b 



Figure 2: Unpolarized decay rate for model 1 as a function of y = m H +/mt (Fig. 2a; tan/3 = 10) 
and as a function of tan/3 (Fig. 2b; m H + = 120 GeV). Other parameters are rrib = 4.8GeV, 
m t = 175 GeV. The barely visible dotted lines show the corresponding mb —> 0 curves. 


In model 1 the y dependence of the unpolarized and polarized rate is essentially de¬ 
termined by the (1 — y 2 ) 2 dependence of the overall factor T in Eqs. flQali and flOblh The 
overall scales of both the unpolarized and polarized rates are set by cot 2 /3, which means 
that their ratio, the polarization asymmetry «#, is independent of the value of tan/3. 

In Fig. 2a we show the unpolarized rate as a function of y — mH+/m t for tan/3 = 10, 
which exhibits the (1 — y 2 ) 2 dependence of T in Eq. (fTOl) . Compared to the partial Born 
term rate r i ^ l y++fe = 1.56 GeV the rate into a charged Higgs is generally quite small 
except for small tan/3 values. This can also be seen in Fig. 2b where we plot the rate 



Figure 3: Polarization asymmetry an for model 1 with mb = 4.8 GeV and mt = 175 GeV as a 
function of mjj+/mt- No value of tan/3 is given since an does not depend on tan/3 in model 1. 
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Figure 4: Unpolarized decay rate for model 2 as a function of y = m>u+/mt (Fig. 4a; tan (3 = 10) 
and as a function of tan/3 (Fig. 4b; m^+ = 120GeV). Other parameters are mb = 4.8GeV, 
mt = 175 GeV. The barely visible dotted lines show the corresponding mi, —► 0 curves. 


as a function of tan/3 for a sample value of m#+ = 120 GeV, which shows the 1/ tan 2 /3 
fall-off behavior of the overall factor (a 2 + b 2 ) oc 1/tail 2 /? in the rate expression Eq. ([9ajh 
The rate into a charged Higgs can be seen to be quite small for tan/3 values exceeding 
tan/3 = 2. One finds equality of the partial rates into a W + and H + only at tan/3 = 0.56 
for rn H + = 120 GeV. Such a small tan /3 value is excluded by the indirect limits in the 
( rriH± , tan /3) plane [16] . 

In Fig. 3 we display the behavior of the polarization asymmetry an as a function 
of y = mn+jmt. As expected from the flat behavior of (a 2 + b 2 )/{2ab) the polarization 
asymmetry an stays very close to its maximal value au ~ 1 except for the endpoint 
region y —> ( 1 — e). At y — 1 — e the polarization asymmetry drops to a# = 0 (not 
shown). We do not show a plot of the polarization asymmetry as a function of tan/3, 
because, as mentioned before, the polarization asymmetry does not depend on tan/3 for 
model 1. 

The y and tan /3 dependence of model 2 is somewhat more involved. The overall scales 
for the unpolarized and polarized rates are now set by (cot 2 /3 ± e 2 tan 2 /3), respectively. 
The y dependence of both the unpolarized and polarized rates is again close to the e —> 0 
limit since they are dominated by the (1 — y 2 ) 2 behavior of T. However, in the case 
of the unpolarized rate, the y dependence is slightly modified by the curly brackets in 
Eq. (l9a]h In Fig. 4a we display the dependence of the rate on y = mn+ /m t where we take 
again tan /3 = 10 as a numerical example. As in model 1 the functional behavior is again 
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Figure 5: Polarization asymmetry an for model 2 with mb = 4.8 GeV and mt = 175 GeV as a 
function of m H +/mt (Fig. 5a, tan/3 = 10) and as function of tan/3 (Fig. 5b, m H + = 120GeV). 
The dotted lines (barely visible in Fig. 5b) show the corresponding m& —> 0 curves. 


essentially determined by the overall (1 — y 2 ) 2 behavior of Eq. ([9a]) except for the region 
close to the endpoint where the second term in the curly brackets of Eq. (l9a]j comes into 
play. The rate is largest for m H + = 0 and drops to zero towards the phase space boundary 
y — 1 — e. The model 2 rate is considerably larger than the model 1 rate, depending, of 
course on the value of tan/3. This is because of the overall rate factor (a 2 + b 2 ) whose 
model 2/model 1 ratio 

(a 2 + b 2 ) | model 2 (1 + e 2 tan 4 /3) 

(a 2 + & 2 )| mod eii ~ (1 + e 2 ) 1 j 

is larger than 1 for tan/3 > 1. For example, for tan/3 = 10 and e = 0.027 one has 
(a 2 + & 2 ) |model 2 /{° 2 + b 2 ) | m odei i = 8.28. The second term in the curly brackets of Eq. (I9al) 
gives a much smaller enhancement. For example, for y — 0, tan/3 = 10 and e = 0.027, 
the enhancement due to the second term in the curly brackets of Eq. (I9al) is only 3.5%. 
Nevertheless, the model 2 rate is still small compared to the Standard Model Born term 
rate T t ^ w + +b — 1.56 GeV for most of the range of tan/3 value. 

In Fig. 4b we plot the tan/3 dependence of the unpolarized rate for m#+ = 120 GeV. 
The tan /3 dependence of the rate is determined by the overall tan /3-dependent factor 
(a 2 + b 2 ) oc (1 + e 2 tan 4 /3)/tan 2 /3 with a slight modification due to the second term in 
the curly brackets in Eq. (jUa|) . As in model 1 the rates drop very fast at the origin for 
small tan/3 values due to the dominant 1/ tan 2 /3 behavior of the overall factor (a 2 + b 2 ). 
Contrary to model 1 the rate rises quadratically for larger tan /3 values after going through 
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a minimum at tan/3 ~ 7 due to the tan 2 (3 term in the overall rate factor. It reaches 
equality with the Standard Model Born term rate Y t ^ w + +b = 1.56 GeV at tan/3 = 0.56 
and tan/3 = 64.83. Again these values of tan/3 are excluded by the indirect limits in the 
(m#±, tan /3) plane [16] • 

All of what has been stated in the last paragraph does not change qualitatively but 
quantitatively if one uses a running mass for the bottom quark instead of a fixed pole mass. 
At one-loop running one would then have m b {m t ) = 1.79 GeV and e = rh b (m t )/m t = 0.010 
instead of e = 0.027. The main effect would be a reduction of the model 2 rate for larger 
values of tan/3 due to the overall factor (a 2 + b 2 ). For example, for tan/3 = 10 one would 
have a rate reduction by a factor of 0.24 due to the factor (a 2 + b 2 ) when one uses a 
running bottom quark mass instead of a fixed pole mass. 

The main qualitative features of the behavior of the polarization asymmetry an in 
model 2 can again be understood from the m b —> 0 behavior of an given by the ratio of 
Eq. (l9bl) and Eq. (l9ajh One has 

1 — e 2 tan 4 /3 f 4e etan 2 /3 1 1 
m™o aH = 1 + e 2 tan 4 /3 { + (1 - y 2 ) l + e 2 tan 4 /3j ' ^ ^ 

The y dependence comes only from the second term in the curly brackets of Eq. (fT2l) and 
is therefore not very pronounced. This can be seen in Fig. 5a, where we show a plot of 
the polarization asymmetry an as a function of y for tan/3 = 10 in model 2. an is large 
and negative with only little dependence on the Higgs mass except for the region close to 
the phase space boundary y — 1—e where an approaches zero, an is negative due to the 
fact that the numerator factor (1 — e 2 tan 4 /3) is negative for tan/3 = 10. For the smaller 
y values the functional behavior is dominated by the first term in the curly brackets of 
Eq. (fT2]h and therefore the curve is rather flat. Towards y — 1 — e the second term in the 
curly bracket of Eq. (1121) becomes dominant and the curve bends up and finally reaches 
zero at y — 1 — e. The e = 0 approximation can be seen to be quite good up to y ~ 0.7. 
In Fig. 5b we show a plot of the polarization asymmetry an as a function of tan/3 for 
m H + = 120 GeV, i.e. y = 0.69. Contrary to model 1, an shows a strong dependence 
on the value of tan/3, an is positive/negative for small/large values of tan/3 and goes 
through zero for tan /3 = \fm^Jm b = 6.04 close to where the rate has a minimum. Beyond 
the zero position an rapidly approaches values close to an = —1- Both Figs. 5a and 5b 
show that m b —> 0 is a very good approximation. At the scale of the figures the dotted 
curve cannot be discerned from the exact full curves. 
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3 Virtual corrections 


The virtual one-loop corrections to the (tH + b)-ve rtex exhibit both IR- and UV-singularities. 
As mentioned before, the UV-singularities are regularized in D = 4 — 2 cn dimensions, 
whereas the IR-singnlarities are regularized by a small gluon mass m g . The scaled gluon 
mass is denoted by 



(13) 


m t 


The renormalization of the UV-singularities is done in the “on-shell” scheme. 

The renormalized amplitude of the virtual corrections in the right- and left-chiral 
representation can be written as 


Mioop = u b | (al + 675 ) Ai + aA 2 + SA^u t , 


where the functions Ai and A 2 read (Cp = |) 


(14) 


Ola 


Ai = ^-C F { (ml + ml - mjj) C 0 - 2rnf - m 2 H + m b (m t + m b ) 

Z7T 


2ml - m 2 H + m t (m t + m b ) C 2 + 2 B 0 — 1 >, 


C 1 


a* 


A 2 = — C F 
' 2tt 


$y2m b m t {C 1 + C 2 ) 


(15a) 


The standard one-loop integrals B 0 , C 0 , C\ and U 2 are given in Appendix [AI After a 
linear transformation these loop functions can be seen to agree with Eqs. ( 8 ) and (9) in 
Ref. [Tj. Following Refs. 0 0 US E! the counter term of the vertex is given by 


iA = (“+^(^-1) + ^-!)-^) 

+ no 

In the “on-shell” scheme the wavefunction renormalization constant Z\ and the mass 
renormalization constant 8m q can be calculated from the renormalized QCD self-energy 
S 9 (p) of the quarks q — t,b. The evaluation of the two conditions T, q (p)\^ =rriq = 0 and 
dT,q(p)/d$\tf =mq = 0 lead to the following renormalization constants 
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z l ,\ = 

= 1-^C F 


An 

5m a 

= —0 Fin, 

H 

An 


-7 E + ln 


47TyU 2 


m: 


UJ 


m 


-2 In —^-+4 


q ™q 

2 


3 47T/T 

' Q1 " +4 , 


a; 


mn~ 


(17a) 

(17b) 


where p is the four-momentum and m g is the mass of the relevant quark. 

Putting everything together the virtual one-loop contributions to the unpolarized and 
polarized rates read (A = A(1 ,y 2 ,e 2 )) 


r i 0 o P = r Born (zl - i + z b 2 - i - — - — + 2A, 

V mt m b 

+ ^ (2 a 2 A 2 ( ( i + e ) 2 - y 2 ) - 2a6 (i - j, 2 + < 2 ) ( 

rL_ = r£..fzj-i + z*-i -—+ 2A, + a 2 

m b 


5m t 5m b 
nn t m b 


loop 1 Born l ^2 




167T 


V m b m t 


(18a) 


The renormalized virtual one-loop correction to the unpolarized rate f!18ap is in agree¬ 
ment with Ref. [7]. The result for the virtual one-loop correction to the polarized rate is 
new. Note that the infrared divergent terms residing in the renormalization factor Z\ and 
in the integral term C 0 in are proportional to the Born term rates T Born and T^ orn , 
respectively. 


4 Tree—graph contributions 

The 0(a s ) real gluon emission (tree-graph) amplitude reads 

2 p° - ^ 7 ® 2 pi + 7 ^ 


. , _ A® _ 

J^Ltree 9s ^ ^b 


al + 675 )u t £l(k, A), 


(19) 


2k-p t 2 k-p b 

where the first and second terms in the curly brackets refer to real gluon emission from 
the top quark and the bottom quark, respectively. The four-momenta of the charged 
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Higgs and the gluon are denoted by q and k. The polarization vector of the gluon with 
momentum k and spin A is denoted by e(k, A). 

In the square of the tree-graph amplitude the terms without gluon momentum in 
the numerator, i.e. the terms proportional to pt and pb in the numerator lead to IR- 
divergences, which are regulated by a small gluon mass. Separating these terms will 
split, the squared tree-graph amplitude into an IR-convergent part and an IR- 

divergent part proportional to the soft gluon factor \M.\ 2 SG f as follows: 

\M tree \ 2 = \M?™\ 2 + \Mq\ 2 \M\ 2 sgf (20) 

where the factor |.Ado | 2 and the universal soft gluon (or eikonal) factor \M.\ 2 SG f are gi yen 
by 


1-^lsGF — 

\M 0 \ 2 = 


—A-KasCp 


m: 


+ 


mr 


- 2 - 


Pb-Pt 


_(k-p t ) 2 ( k-p b ) 2 (k-p b )(k-p t ) _ 

2 (p t ■ Pb)(a 2 + b 2 ) + 2(a 2 - b> 2 )m t m b + 4 abm t (p b ■ s t ). 


( 21 ) 

( 22 ) 


Note that the factor |2Vdo | 2 has the same analytical form as the squared Born term 
amplitude |2V4o | 2 listed in Eq. ([ 6 ]) but must of course be evaluated for p t = pb + q + k, i.e. 
one has \Mo\ 2 (k = 0 ) = |-M 0 | 2 . 

When integrating \M.o\ 2 \J\A\ 2 SGF one has to take care of the fact that both terms in the 
product depend on the gluon momentum. Since (|WTo1 2 — |-Mo| 2 )|A4|| GF is convergent and 
thus can be integrated without a gluon mass regulator, we further isolate the IR.-divergent 
part by writing 


I M tree \ 2 = {\M%Z I 2 + (lM.1 2 - |Vfol 2 )|X|| GF } + |A<ol 2 |Vf|| G F. (23) 


where the gluon momentum dependence in the IR divergent last term of (ITU) solely resides 
in the soft gluon factor |A4|| GF . 

The same universal soft gluon factor |A4|| Gi? appears in the calculation of the radiative 
corrections to t —> W + + b. We can therefore take the result of its phase space integration 
(with m g 7 ^ 0) from Eq. (63) in Ref. [14] , 

The IR-convergent part of \A4 tree \- is thus given by 
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| M 


conv 1 2 
tree I 


+ \M 0 \ 2 - \Mq\ 2 )\M\s GF = 8 tt a s C F 


k-q 


{k-p t )(k-p b ) 


(a 2 + 6 2 )(fc- q) + 


-2 abnrii 


m: 


mi+m 1 H+ — mi 


m 


H+ 


2k-p t 


mi 


m 


H+ 


k-pt k-p b 


k-q 


k-q 

(q-s t )~(q-s t )\k=o 


(k-s t ) + (q-s t ) 


( 24 ) 


It is quite remarkable that the unpolarized and polarized pieces of the convergent 
tree-level contribution are proportional to (a 2 + b 2 ) and 2 ab, respectively. 

The phase space integration is done with respect to the gluon energy ko and the H + 
boson energy E H + 1 where the k 0 integration is done first. It is convenient to introduce the 
scaled invariant mass of the bottom quark-gluon system z = 1 . Using the relation 

E h + = (m 2 + m 2 H+ — m 2 z)/{2m t ) we change to the integration variable z. 

Note that the last line of the Eq. (124)) comes from the term (jA^ol 2 — |A4 0 | 2 ) |A4|| GF 

in Eq. (1231) . The relevant scalar products are evaluated in the top quark rest frame. They 
read 


<«■*> =-5 
(q-s t )\k=o = — 2 


a/a(T , y 2 , z ) m t cos 9 P , 
\/A(l,y 2 ,£ 2 ) m t cos 6 P . 


(25) 

(26) 


After ko integration the remaining phase space integrations can be reduced to the 
following classes of basic integrals (A' = A(1 ,y 2 ,z)): 


R(n) 



R(m, n ) 



S(n) := 


S(m, n) 


[ dz - 1 —= 

J (z — e 2 ) a/A" 

f z m 

:= / dz _ In 

./ [ 


\l-y 2 +z + VX' 

In - -= , 

.1— y 2j rz— V^V. 

1 - y 2 + z + V*- 
1 — y 2 + z — v^V. 
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(27a) 

(27b) 

(27c) 

(27d) 






















with the integration limits e 2 ^ z ^ (1 — y) 2 . The basic integrals needed in the present 
application are listed in Appendix O 

Finally, the nnpolarized and polarized tree-graph contributions read 


p 

L tree 


Anmt 


^C F m 2 (a 2 + b 2 ) | \r(0 , -1) + Ai_A> -1) 


1 y2 + 3t2 m-l, -1) + i« 2 S(0,0) - is(l, 0) !■ + PS 2 -‘ r B „„ 5(A) 


(28a) 


and 


r = 

tree 


Airm t 
1 


^ C F m 2 2 ab \ - 2^X R(- 1) + 2A R( 0) + -(1 - y 2 fe 2 R{-2 , 0) 


An 


j ((1 - V 2 ) 2 + 2 (3 - y 2 ) e 2 ) 0) - 1(2 + 10 y 2 - e 2 ) K(0,0) + Ifl(l, 0) 

+(1 - V 2 +1 2 ) v^SIO) - (1 - !, 2 + e 2 )A 5(1) + \/A5(0,0) + i(4y 2 (1 - j, 2 ) 


e 2 - 2 e 4 


(7 + 5 y 2 ) 

2 x Born S(A) 


) 5(0,1) - 1 (3 - 3 y 2 + e 2 ) 5(1,1) - |s(2,1) 


+pSo 1 r p 


(28b) 


As before, PS 2 = \/A/167r m t is the two-body phase space factor. 5(A) is the integration 
of the soft gluon factor |Ad|| GF defined in Eq. (I2TT) . i.e. 


-S'(A) := (-47rm t ) 


1 1 

2 m t (27r) 5 


d 3 p b d 3 q d 3 k 
2 Eb 2 Eh 2Eg 


^ A \pt ~ Q ~ Pb — k)\AA \ 2 sgf . 


The integration 5(A) was done e.g. in p2E]. For completeness we list the result in 
Appendix [B] 
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5 Results 


For ease of comparison with the results of 13 E] we define the following abbreviations: 


p 0 = -(l + e 2 - 2 / 2 ), 

Ps = \\/K 1 p 2 ^v 2 ), 

p±=P0±p3, 

V 1 I P+ 

Y P = o ln — > 

2 p_ 


= ^(1 — e 2 + 2 / 2 ), 

w 3 =P3, 
w± = Wo ± u>3 , 

V' _ 1 1 

Yw o m , j 
2 w _ 


(29) 


where the hat symbols on the scaled momenta pi and tty accentuate the fact that we 
are dealing with dimensionless quantities. The complete 0(a s ) results are obtained by 
summing the Born term, the virtual one-loop and the tree-graph contributions, i.e. 


= hr+ pr p cosOp) 

a cosOp 2 


Born 


+ r 


NLO 


+ (r F Born + r 


NLO 


P cos 6 f 


The unpolarized 0(a s ) corrections are given b}0 


T 


NLO 


r loop t Tj ree 


a s 

8n 2 


C F m t 


(i a 2 +b 2 ) G + + (a 2 — b 2 ) eG_ + abG 0 


with the coefficient functions 


r9 


API 


G + =p 0 n + Pops ^ - 4ln(d^) + — Yp(2 - y 2 - 4 y A 
L2 ey J 4 y z 

+3 y e - 2e 2 - 2e 4 + 2e 6 - 4 y 2 e 2 - 5 y 2 e A ), 

G_ =W +P3 
Go = - 6 p 0 p 3 lne, 


6-4 ln(—)] + - y 2 - 2e 2 + e 4 - 3y 2 e 2 ), 

ey J r 


where 


n = Apo 


p-s , ,, ,._/ 4 2/P3- 


Li 2 (p+) - Li 2 (p_) - 2Li 2 (l - —) + Tpln( 

p + p% 

+21^(1 — e 2 ) H— 2 ^ 3(1 + y 2 ~ A) lne . 

y 


- Y w ln e 


(30) 


(31) 


(32) 


( 33 ) 


1 This result is written in the same form as in [6|. 
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The polarized 0(a s ) corrections are given by 


Y'Pol _ t~iP _i_ rP 

1 7 VXO — 1 i 1 / 


loop ' tree 




87T 2 


m t C F l — 3(gT + lne + ab 


11 + 28 y- 16 y 2 - 8y 3 + 7 y 4 


+e 2 (4 + 8 y- Uy 2 ) + 7e 4 ) + (2 - 9 y 2 + y 4 - e 2 (4 + 3y 2 ) + 2e 4 ) |§ 

+ 8 ^ ln ( ( 1 _ 1 y )2 y _, e2 ) + ( 3 - 3 P + 2e2 ( 4 + y 2 ) ~ 2e4 ) ln (“7^) 


1 -v 

n l(T 

8 Pi™o /+t: 1 — 2/s OT . 1 — 2/s 




r 


lne + 4p 0 p 3 2Li 2 (l- 7 —) - 2Li 2 (l - 

p_ 


P+ 


(1 — V 2 ) — e 2 

-Li 2 ( , u>_) + Li 2 (++) + 21n(-- - - )Y p 


-(2 +y 4 - e 2 (3 + 2 y 2 ) + e 4 ) (2Li 2 (y) - Li, 2 (uL) - Li 2 («>+)) 


(34) 


As mentioned before the 1 R and mass singularities cancel in the sum of the one-loop 
and the tree-graph contributions. 

Next we discuss various limiting cases for the unpolarized and polarized rates, which, 
among others, serve to check on the correctness of our results. We have checked that our 
TTib 7 ^ 0 results for the unpolarized rate agree with those given in [ 6 ]. We do not, however, 
agree with the m b 7 ^ 0 results of [7]. 

The limit m F + —> 0 is of interest since, according to the Goldstone equivalence theo¬ 
rem, the unpolarized and polarized rates for t —> H + + b become related to the unpolarized 
and polarized longitudinal rates of t — > W + + b in the limit m F + ,rnw+ —► 0. For model 1 
and tan fi — 1 in the m F + —> 0 limit one has 


rri'f Gp 


ot s 


7 r 


lim r = -l^ = \V tb \ 2 (l-ey<l + —C 

m H+ ^0 87T y/2 tb ' ' 

—2Li 2 (l — e 2 ) — 2 


2 1 — e 2 (1 — e 2 

0 - 77 hi 

1 + e 2 


1 + e 2 [5 - 22e 2 + 5e 4 
1 — e 2 L 4(1 — e 4 ) 

4 — 5e 2 + 7e 4 


2 lneln(l - e 2 ) 


(l-e 2 )(l-e^) 


ln e 


(35a) 


lim r p = m ‘ Gf 

m H+ _> 0 


87T y/2 


a. 


V tb \ 2 (l-e 2 y{l + -C F -— 2 

7 r 1 — G 


- ^(5 + e 2 ) — 2(1 + e 2 ) lneln(l — e 2 ) 


- 2(1 - e 2 )In 


1 — e 5 


4 + 5e 2 


lne + (1 — 2e 2 )Li 2 (l — e 2 ) 


(35b) 
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We have checked that expressions Eqs. (I35aj) and (I35bh agree exactly with the m\y+ —> 0 
limit of the corresponding longitudinal and polarized longitudinal rates in the process 
t —> W + + b listed in |14j . This is nothing but the statement of the Goldstone equivalence 
theorem. Our unpolarized result in Eq. (I35a.p agrees with the corresponding m#+ —> 0 
result in [b]. 

For the sake of comparison with results in the literature we take both mp+ —> 0 and 
nib —> 0 in Eqs. (T3TT) and (134)) without the cautionary proviso of keeping the term nib tan /3 
in model 2. In both models 1 and 2 one then obtains 


lim T = 

m H +^ 0 


ml Gi 
8n y/2 


\V tb \ 2 cot 2 f3 


1+^C F (--— 
2tt V2 3 


lim V F =^^\V th \ 2 cot 2 (5 


m H +^ 0 


87T y/2 


1 -=!LC F [^-^ 
2 tt V 2 3 


(36a) 


(36b) 


which, when setting cot (3 = 1, agree exactly with Eqs. (48) and (49) of [TT] and the 
e —> 0 limits of Eqs. (I35al) and fl35bl) . The unpolarized rate in this limit agrees with the 
corresponding results in [Slt6[lBj- 

When mu+ approaches m* for m b —■ > 0 one has 


lim —-= 1 + —Cp 

rn H+—>’ rn t I Born 27T 


13 4vr 2 


3 ln(l — y 2 ) 


T p a 

lim —p— = 1 + t ~^C F [1 - 7T 2 - 3 ln(l - y z 

m H +^mt 1 Born 2vr 


(37a) 

(37b) 


where the tan f3 dependence has disappeared by dividing out the relevant Born terms. For 
the unpolarized rate the limiting expression agrees with the corresponding limit given in 

K'l 

Finally, we consider the limit m b —* 0 keeping the charged Higgs mass finite. This 
results in very compact expressions for the unpolarized and polarized rates. Due to the 
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smallness of the bottom quark mass and the fact that the bottom quark mass corrections 
are of 0(m 2 /m 2 ) the m b —> 0 formulae give quite good approximations to the exact 
formulae for Higgs masses as long as the Higgs mass is not close to the top quark mass. 
One obtains 


lim T = 

m;,—> 0 


m t ,2\2(„2 , i2' I - a 


167T 


(1 - y‘)\a‘ + 6 2 )p + 

\ a z + b z 1 — y z 2 tt 


2vr 2 


4 y 1 


i -y‘‘ 


hi y 


r 


-4 In y) ln(l - y ) - 4Li 2 (y ) + 


bf 


+ b 2 


3 In e 


(38a) 


lim T J 

mt ,—>0 


m t 
167r 
2 


(i-yrzab 


1 + £ Cf 


11 — 6 y — 7 y 2 1 + 2 y 2 ^r 


+ 


2(1 + y) 2 (1 — y 2 ) 2 3 


9 y 2 + y 4 


(1 _ y2^y2 

8 + % 4 


ln(l + 2 /) + 


2 — 5j/ 2 


ln(l - y) -AU 2 (y) 


,212 


Li 2 (-2/) 


6 )^ 


2 a 6 


■3 In e 


(38b) 


(i -r 

Again we have checked that the —> 0 limit of the unpolarized rate agrees with the 

corresponding result of [6j|. Note that the seemingly mass singular terms proportional 
to lne in Eq. (I38a,p and (138b|) are not in fact mass singular, since they are multiplied 
by the factor (a — b) 2 , which is proportional to m 2 in both models 1 and 2. Although 
the contributions proportional to (e 2 lne) formally vanish for rri}, 0, they can become 
numerically quite large for rrib = 4.8 GeV in model 2, depending, of course, on the value 
of tan / 3. This can be seen by calculating the ratios of the coupling factor expressions that 
multiply the lne term in Eqs. f!38al) and (138b|) . In model 1 one has (a — 6) 2 /(a 2 + b 2 ) = 
2 e 2 /(1+e 2 ) and ( a—b ) 2 /(2 ab) = 2e 2 /(l—e 2 ), whereas in model 2 one has (a—b) 2 /(a 2 +b 2 ) = 
2e 2 tan 4 /?/(l + e 2 tan 4 /3) and (a — b) 2 /(2ab) = 2e 2 tan 4 /?/(l — e 2 tan 4 /3). In model 1 the 
lne contribution is negligible but not in model 2. For example, in model 2 one finds 
(a — b) 2 /(a 2 + b 2 ) = 1.77 and (a — b) 2 /(2ab) = —2.31 for tan/3 = 10. With a little bit of 
algebra one finds that the NLO corrections in model 2 are in fact dominated by the (e 2 In e) 
contributions for larger values of tan/3 as also noted in [6]. This is evident in Figs. 4a 
and 4b where the radiative corrections to the Born term rate can be seen to be as large 
as —50% compared to the « —10% expected from the corresponding corrections in the 
decay t —> W + + b p3]. Note, though, that the radiative corrections become much smaller 


if one uses the running bottom quark mass mj(mt) = 1.79 GeV, i.e. e 2 lne = —0.00048 
instead of the fixed pole mass nib = 4.8 GeV, i.e. e 2 lne = —0.00271. 

We now turn to a more detailed discussion of our numerical results on the radiative 
corrections to the unpolarized rates and the polarization asymmetry in model 1 and 2. 
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We start our discussion with model 1. As input values for our numerical evaluation 
we use 7rtb = 4.8 GeV and m t = 175 GeV as in the case of Born term. The strong coupling 
constant is evolved from a s (Mz) = 0.1175 to a s (m t ) = 0.1070 using two-loop running. 
Fig. 2a shows that the radiative corrections lower the model 1 rate by approximately 
10% over the whole range of y = rrifj+ /m t . The radiative corrections in Fig. 2b are of 
similar size but cannot be discerned at the scale of the figure. The radiative corrections 
to the polarized rate are of similar size and go in the same direction, which means that 
the radiative corrections to the polarization asymmetry are quite small. This can be seen 
in Fig. 3, which shows that the radiative corrections to the polarization asymmetry are 
indeed quite small and lower an only by ~ 2%o. No value of tan/3 is given in Fig. 3, since 
in model 1 the polarization asymmetry does not depend on tan/3 at LO and NLO, and, 
in fact, at any order in a s . 

Figs. 4a and 4b show that in model 2 the radiative corrections are substantial, which 
is mostly due to the e 2 In e contribution discussed above. They reduce the LO rate by 
~ 50% over much of the shown y and tan/3 ranges. The barely visible dotted curves in 
Figs. 4a and 4b are drawn using the “kinematical” mb —■► 0 approximations Eqs. (j9al) and 
(|9b]i (LO) and Eqs. (I38af and (I38bl) (NLO). In these equations the bottom quark mass 
has been set to zero whenever the scale of rrib is set by rri t as in the kinematical factors and 
not by mt cot 2 /3 as in the coupling factors. As Fig. 4 shows, the “kinematical” mb —> 0 
approximation is an excellent approximation for both the unpolarized and polarized rate. 

In Fig. 5 we show that in model 2 the radiative corrections to the LO Born term result 
are substantial and reduce the size of the polarization asymmetry by ~ 25% over much 
of the range of the Higgs mass. The mb —■► 0 approximation is quite good except in the 
endpoint region. At the scale of Fig. 5 the rrib % 0 corrections are barely visible. In 
Fig. 5b we fix the mass of the charged Higgs boson at m H + = 120 GeV and vary tan/3 
between 0 and 30. The LO zero position tan/3 = \Jm t fmb = 6.04 is shifted upward by 
approximately one unit by the radiative corrections. The radiative corrections can be seen 
to become quite small for the larger tan/3 values. The radiative corrections are largest 
around the zero position of tan /3 at tan /3 ~ 7. 

6 Concluding remarks 

We have calculated the 0(a s ) radiative corrections to polarized top quark decay into a 
charged Higgs and a bottom quark in two variants of the two-Higgs-doublet model. We 
have checked our unpolarized results against known results and found agreement. Using 
the same techniques we have calculated the polarized rate. Further, we have compared 
our polarized results with the corresponding polarized results in the decay t —> W + + b 
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appealing to the Goldstone equivalence theorem. Because of our numerous cross-checks 
we are quite confident that our new results on the polarized rates are correct. We have 
found very compact 0(a s ) expressions for the unpolarized and polarized rates in = 0 
limit, which can be usefully employed to scan the predictions of the 2HDM (mjj+, tan/3) 
parameter space. 

We have found that a measurement of the polarization asymmetry in the decay of 
a polarized top quark into a charged Higgs and a bottom quark t(|) —► b + H + can 
discriminate between the two variants of the two-Higgs-doublet model investigated in this 
paper. In model 1 the polarization asymmetry does not depend on tan (3 and stays very 
close to an = 1 for a large range of charged Higgs mass values. This is different in 
model 2, in which the polarization asymmetry strongly depends on the values of tan (3 
and the charged Higgs mass and can vary between an = 1 and an = — 1 depending 
on the parameter values. The radiative corrections to the polarization asymmetry are 
quite small (« 2%o) in model 1. Again this is different in model 2 where the radiative 
corrections to the polarization asymmetry are important and can become as large as 25%. 
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A Loop integrals 


A detailed discussion of the one-loop integrals can be found in e.g. Sec. 4 of 
The two-point one-loop scalar integrao 


Bo(q, m t , m b ) : = 


A 


4 -D 


d U k 


in 2 J (2 tt) d 4 [(p f - kf - mf\ [( p b - kf - m 2 ] 


= -7 E 

u 


2 +In 


Att/j , 2 


m; 


+ 


y 


e“. 2p 3 


r 


In e -I-- Y p . 

r 


(39) 


2 Note that q = pt — Pb, q 2 = nkj {+ and mt > mb + mg+ . 
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The three-point one-loop scalar integral: 


C 0 (q,m t ,m t ,m b ,m g ) := — 


4 -D 


d U k 


in 2 J (2tt) d 4 [Q 0h - kf - ml] [(p t - k) 2 - ml] (k 2 - m 2 g ) 

— 2 { (lne + Y p ) (Y p + 2 Y w ) + (In e — 2 In A) Y p — Li 2 (l ——) + Li 2 (l — ^—— / 

mf2p 3 { w + P+w+ 


(40) 


The three-point one-loop vector integral: 


k 


C^(q,m t ,m t ,m b ,m g ) := — 


4 " D f d D k 


E L 


in 2 J (27r) D 4 [(p 6 - A:) 2 - m 2 ] [(p t - k) 2 - ml] (k 2 - m 2 ) 

= C\ pi T C'2 p b ■ (41) 


where 


C\ = 


1 1 
mly 2 


Y, 




Co — 


, , lne+(l +y 2 -e 2 )^- 
m% y [ 2p 3 


1 1 
T2T2 


(42) 

(43) 


B Integration of the soft gluon factor 


S{ A) := (-47rm t ) 


Ilf d 3 p b d 3 q d 3 k 

2 m t (27r) 5 J 2 E b 2E H 2 E g 


5 {A \pt - q-pb~ k)\M\ 


SGF 




- 4^ I 4 ^3 


ln( 


4p|. 

67/A 


- 2 


— 2(1 — e 2 )Y w + 2e^ Y p + 2p 0 {2Li 2 (l ——) 

l p + 


T • (1 w- . 2p 3 , 

—Li 2 (1 ——) + Li 2 (-—— 

w + p + w+ 


-Yv 


2 ln( 




P+U ! + A 


-n + 1 


(44) 


where the soft gluon (or eikonal) factor A4 | | CF is dehned in Eq. 
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C Basic integrals for the tree—graph phase space in¬ 
tegrations 

Details of the calculation of basic integrals can be found in [13]. Note that our notation 
differs from the notation in pH] , 


#(-1) = - 2^3 + 2p 3 in( 2w +^ 3 ) - (l - e 2 + y 2 ) Y w , 

y 2 e 2 


R{ 0 ) = - In 


(i -yf 


+ ln( 


W±s 

' ye 2 ' 


2 L(1 + y) 2 — e 2 . 

i?( 0 , 0 ) = (1 - 2/) 2 - e 2 , 

R( o, -i) = P3 (i - e 2 + y 2 ) - Zy 2 Y w , 
fi(-l,0) = 2 Infill), 

«(-l, -1) = -2p 3 + 2 (1 - y 2 ) F p - 2y 2 Y w , 
7?(-2,0) = e - 2 - (1 - s/)- 2 , 


-1) = % - -2-j 

e" 1 1-?/- 


(l+l/ 2 ) + Z/ 2 ^™ 


^( 1 , 0 ) = 




4 A 
2 ’ 


5(0) = Li 2 (w + ) - Li 2 (w_) - 2Li 2 (l - t) +21n(^±^)F p - 2F p 2 + 21^ lire, 

p+ y£2 


5(1) = -x— { Li 2 (—^—) — Li 

P3 < w + 


yp~ 


p+w+ 


) + in (—)y„-yA, 

yt 2 J 


5(0, -1) = -1 [(1 - V? - f 2 ] [ £ 2 - (3 - y) (1 + y) ] - (1 - y‘) 

+2y 2 2Li 2 (y) - Li 2 (ro_) - Li 2 (w + ) + 2 (l - e 2 + y 2 ) p 3 F p , 

5(0,0) — 2 (p 3 — e 2 Y p — y 2 Y w ) , 

5(0,1) = Li 2 (w_) + Li 2 (w>+) - 2Li 2 (y), 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 
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(59) 


S{ 1, 0) = P j- (1 + e 2 + 5 y 2 ) - e% - y 2 ( 2 + y 2 ) Y w , 

S(l, 1 ) = e 2 - (1 - y) 2 + 2 (1 - y 2 ) ln(^—-) - 4 p 3 Y p 

+ (l + y 2 ) Li 2 (w_) + Li 2 (D+) - 2Li 2 (y) , (60) 

S(2,1) = j [e 2 - (1 - J,) 2 ] [4 + £ 2 + (1 - yf + 8 S 2 ] + 3 (l - y‘) ln(i^) 

- (l + 4 y 2 + y 4 ) 2Li 2 (y) - Li 2 (w_) - Li 2 (w + ) 

-2p 3 (3 + e 2 + 3y 2 )Y p . (61) 
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